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(-H I Abstract 

We use WKB Approximation in the calculation of slope and curvature of Isgur-Wise 
function in a QCD inspired potential model. This work is an extension of the approximation 
r~| I methods to the QCD inspired potential model. The approach hints at an effective range of 

distance for the calculation of slope and curvature of Isgur-Wise function. Comparision is also 
made with those of Dalgarno and VIPT to see the advantages of using WKB Approximation. 
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1 Introduction 

Probably this is the first attempt to use WKB Approximation (WKBA) for the calculation of slope 
and curvature of Isgur-Wise (IW) function with linear cum Coulombic pottential of QCD. The 
WKBA which gives approximate but direct solution of Schrodinger equation is applicable to re- 
', gions where potential energy variation is small. There are few works [U El |3], where linear cum 
Coulomb potential is used under WKBA to calculate the energy spectrum , decay width etc. With 
the success of those models ,we here try to extend WKBA in the calculation of slope and curvature 
of IW function to sort out the advantages as well as the limitations of the method. 
It is observed that so far the results for slope and curvature are not satisfactory for Be meson [4- 
10]in our QCD inspired potential model with significant confinement effect b = O.lSSGeV"^ [Hj.We 
opt this method in the hope of doing better for the meson. 
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2 Formalism 

The one dimensional (say x) Schrodinger equation is [19] : 

fl + 2f,[E-V{x)]^ = 



For E > V{x) ,it becomes : 

f^ + e{x)^ = o (2) 

where 

k^x) = 2fi[E- V{x)] (3) 

The WKB solution of Eq.(2) is : 



-.Cos{ J k{x)dx - ^) (4) 



'k(x 

For E < V{x) , the WKB solution is either: 

^(a;) = ^=e- (5) 
./t{x) 

or 

^{x) = ^!=e/^' (6) 

where 

k'\x) = 2fi[V{x)-E] (7) 
and A^,A^' , A^" are the normalization constants . 

The WKB quantization condition for the allowed bound state energy is [19]: 

(^rir + ^) = / k{x)dx (8) 
where A and B represent the turning points at which E = V{x) and n,. > 0. 

For the three dimensional radial Schrodinger equation which contain an additional centrifugal 
potential to the original potential y(r),the independent variable r is converted to x to give 
the three dimensional equation a one dimensional form and in the process we instead of the integral 

k{x)dx are left with the integral given as [19] : 



k(r)dr = [E- V(r) - ^ ^ dr (9) 
A JA \ 2fir^ J 

We note that the term -^j^ has changed to ^^^2 and the WKB solution is applicable to radial 
Schrodinger equation like the one dimensional Schrodinger equation in the variable x. 



With the original linear cum Coulombic potential [12] the WKB quantization condition in this 
case (Eq.9 above) for the ground state (n^ = 1, / = 0) is : 



1\ /„ , a 1 Y 
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The limits A and B are the turning points given by the positive roots of the cubic equation: 

br^ - Er^ -ar + - = (11) 
4 

We also note that 

a = (12) 

At the turning points, the WKB solutions become invalid and we find wave functions only in 
the regions r < A{E < V),A < r < B{E > V) and r > B{E < V).We use the leading order 
expressions for energy given by Eq.(ll) of Ref[2] with W,^ replaced by E and b respectively: 

The solutions for the respective regions are [19] : 

vl/,(r) = -^e-/'=''^^r<A (14) 

No r rc 

^2(r) = ^ Cos{ kdr - -),A <r < A (15) 
2\/7rkr J 4 



^3(r) = ^ e/^''^^r > B (16) 
2V nk'r 

where Ni,N2 and are the normalization constants in the respective regions while k' and k 
are given respectively by: 



k' = v/2M^ -E) = ^/2^ (br--Er-^ + ^ ) (17) 



and 



k=^2f,iE-V) = \l2f,[Er-br^ + ^-^ 



As we have seen that there are three wave functions corresponding to the three different regions 
so we expect three IW functions also for these regions. 

The IW function and its slope and curvature can be calculated as [20]: 



£,{y) = / 47rr |'i/'(r)| cosprdr (19) 
Jo 

= l-p^y-l) + C{y-lf + ... (20) 

where 

/ = 2/^2 (y - 1) (21) 

We however in this work restrict ourselves to a finite maximum of r (i.e.rmax) instead of taking it 
as 00 in Eq. (19). This is because for infinite maximum of r we have diverging results. Consideration 
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Figure 1: Variation of V{r) vs r for the energy E given by Eq.13.The points of intersection 
between y(r)(sohd curve) and E (dashed curve) are the turning points A and B -the positive 
roots of Eq.ll. 

of finite maximum of r further hints at the usuable range of r in the calculation of IW func- 
tion. However, only for the region r > B, this consideration comes into existence .The regions can 
be located from fig.l. 

The integration involved in different equations (14), (15), (16) lead to the well known elliptic 
integrals [2]. Instead of using the elliptic functions we have integrated them numerically using 
mathematica software. 



We use both the MS [4J and y- schemes [T7 t [T8] to see the role of in WKBA in calculating 
the IW function. 



3 Calculation and Result 

The turning points A,B as well as the Vmax are listed in table 1. Table 2 shows the slope and 
curvature C under MS with and without relativistic effect while table 3 shows the same for V- 
scheme.In table 4 we have recorded the predictions of the other models to make comparision with 
our work. 
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Table 1: List of the turning points A,B in GeV ^. 
Mesons MS scheme V scheme 





A 


B 


A 


B 


D 


0.431 


7.52 


0.337 


7.79 


Ds 


0.353 


7.53 


0.269 


7.8 


B 


0.461 


7.36 


0.444 


7.4 


Bs 


0.374 


7.375 


0.357 


7.41 


Be 


0.195 


7.39 


0.182 


7.43 



Table 2: Values of and C under AfiS-schcme. Subscript 'rel' refers to the inclusion of relativistic 



effect. 


The maximum 


finite distance ta 


ten for all the mesons are 




= 9.5 GeV-^. 






Me- 
son 


Region r < A 


Region A < r < B 




Region 


r > B 






C 


Prel 


Crel 


P' 


C 


Piel 


Crel 


P' 


C 


Prel 


Crel 


D 


3.2x 


4.92 X 


1.2x 


1.61X 


0.357 


4.34X 


0.142 


1.8x 


5.09 


0.065 


5.06 


0.0643 




10-3 


10-6 


10-3 


10-6 




10-3 




10-3 










Ds 


3.78X 


6.66x 


1.51X 


2.33X 


0.533 


0.012 


0.190 


4.3x 


9.07 


0.205 


9.02 


0.2038 




10-3 


10-6 


10-3 


10-6 








10-3 










B 


4.56X 


1.03 X 


2.8x 


5.96X 


0.416 


6.87X 


0.258 4.26X 


6.41 


0.106 


6.38 


0.105 




10-3 


10-^ 


10-3 


10-6 




10-3 




10-3 










Bs 


5.8x 


1.62X 


3.72X 


9.74X 


0.733 


0.0237 


0.436 


0.014 


12.6 


0.405 


12.5 


0.404 




10-3 


10"^ 


10-3 


10-3 


















Be 


1.24x 


6.46 X 


9.1x 


4.4x 


3.07 


0.712 


1.52 


0.353 


90.8 


21.1 


90.5 


21.0 




10-2 


10-^ 


10-3 


10-^ 



















Table 3: Values of and C under V^-scheme. Subscript 'rel' refers to the inclusion of relativistic 



effect. The maximum finite distance taken for all the mesons are r 



9.5 GeV- 



Me- 
son 




Region 


r < A 




Region A<r < B 




Region 


r> B 




P' 


C 


Prel 


Crel 


P' 


C 


Prel 


Crel 


p' 


C 


P'iel 


Crel 


D 


1.87X 


1.62X 


1.66x 


8.54X 


0.296 


3.76X 


0.032 


4.07X 


5.33 


0.07 


5.27 


0.07 




10-3 


10-6 


lo-'i 


io-« 




10-3 




10-^ 










Ds 


2.1x 


1.99X 


2.48X 


1.4x 


0.432 


9.75 X 


0.036 


8.15X 


9.5 


0.215 


9.41 


0.213 




10-3 


10-6 


10-3 


10-^ 




10-3 




10-^ 










B 


4.08X 


8.28X 


2.26X 


4.25X 


0.431 


7.11X 


0.244 


4.0x 


6.46 


0.107 


6.43 


0.106 




10-3 


10-6 


10-3 


10-6 




10-3 




10-3 










Bs 


5.1x 


1.25X 


2.98X 


6.73X 


0.677 


0.022 


0.354 


0.011 


12.7 


0.41 


12.6 


0.407 




10-3 


10-^ 


10-3 


10-6 


















Be 


1.2x 


5.7x 


8.2x 


3.6x 


3.18 


0.74 


1.32 


0.31 


91.5 


21.2 


91.1 


21.2 




10-2 


10-^ 


10-3 


10-5 
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Table 4: Predictions of the slope and curvature of the IW function in various models. 

Model Value of Value of curvature C 



Le Youanc et al 121] > 0.75 



Le Youanc et al [22] > 0.75 > 0.47 

Rosner YM 1.66 2.76 



Mannel [Ml ES] 0.98 0.98 

Pole Ansatz [26] 1.42 2.71 

MIT Bag Model [28] 2.35 3.95 

Simple Quark Model [27] 1 1.11 

Skryme Model [25] 1.3 0.85 

QCD Sum Rule [26] 0.65 0.47 



Relativistic Three Quark Model [24] 1.35 1.75 

Infinite Momentum Frame Quark Model [23] 3.04 6.81 

Neubert [33] 0.82±0.09 

UKQCD [sg 0.83l}f;i 

CLEO [33] 0.76 ±0.16 ±0.08 
BELLE [31 0.69 ±0.14 



4 Discussion and Conclusion 

The method used for the first time in the calculation of IW function results in slope and curvature 
which is comparable with other models for the region A < r < B only.For the other two regions 
r < A and r > B the results are far from being satisfactory. For the region r < A {r > B) the 
values of and C are quite small (large) which shows that the IW function is dependent on the 
distance(r) considered and the slope and curvature increase with r. Further, our assumption of 
finite maximum (rmax) distance instead of infinity as upper limit (r = oo) is also valid because 
with increase in r the slope and curvature of IW function increase to great extent leading to quite 
unsatisfactory result as stated above. Consideration of only the leading term of Eq.ll of Ref[2] may 
be a cause of these results and hence exploration of approximate expression for energy (Eq.l2 of 
Ref[2]) is another prospect in the model. 

We have also noticed that and C decrease with the increase in as as our results are larger 
with MS'-scheme which have smaller than V^-scheme having larger a^. Further, the relativistic 
effect has reduced the values of p^ and C as observed earlier [4- 10]. However, it is interesting to 
observe that for the region r > B, there is a very little effect of as as well as relativistic consider- 
ation introduced through e hinting at the necessity of an effective finite distance as maximum for 
the calculation of IW function . 

To conclude , this work opens the idea of a usable range of distance in the calculation of IW 
function and this can be the alternate consideration that one may take account in renovating the 
earlier work. 
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